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This paper presents revised calculations for the Maki parameters ki and k,2 and the pair potential 
A(r) of s-wave type-II superconductors near the upper critical field H C 2 with arbitrary impurity 
concentration. It is found that Eilenberger's well-known results on «2 [Phys. Rev. 153, 584 (1967)] 
are not correct quantitatively, which are modified appropriately. Calculations are also performed 
for a two-dimensional system with an isotropic Fermi surface. The results in the clean limit differ 
substantially from those in three dimensions with an spherical Fermi surface. This fact indicates 
the necessity of considering detailed Fermi-surface structures for a quantitative understanding of 
the parameters. The coefficient of A(r) oc (H& — B) 1 ^ 2 , which is basic to any theoretical evaluation 
of the thermodynamic and transport properties near is obtained accurately. 



PACS numbers: 74.25.-q, 74.25. Op 



I. INTRODUCTION 

Following the preceding studies; 1 : 2 : 3 : 4 : 5 : 6 : 7 : 8 : 9 : 10 : 11 
Eilenberger— performed an extensive calculation of the 
parameters Ki(T) and K2(T) introduced by Maki^ to 
distinguish temperature dependences of the upper crit- 
ical field H C 2 and the initial slope of the magnetization 
dM/dH, respectively. Based on the s-wave pairing with 
a spherical Fermi surface and taking both s- and p-wave 
impurity scatterings into account, he clarified a basic fea- 
ture that K2 > Ki > ^gl, where kgl is the Ginzburg- 
Landau parameter near T c> i&ii He also found a large 
dependence of the parameters on the p-wave scattering 
strength. This study is undoubtedly one of the basic 
works in the field and has been referred to frequently in 
analyzing experimental results on the quantities. It will 
be shown, however, that his results on K2 are not correct 
quantitatively due to a couple of inappropriate approxi- 
mations adopted. 

This fact also tells us that we are still far from a quan- 
titative description of type-II superconductors. The pa- 
rameter K2 is such a basic quantity that it is relevant to 
all thermodynamic and transport properties near H C 2. 
Indeed, changes of those quantities through H C 2 are pro- 
portional to the spatial average (|A(r)| 2 ) of the pair po- 
tential A(r), and (|A(r)| 2 ) is directly connected with k 2 , 
as seen below. Thus, an absence of a reliable theory on 
K2 also implies no quantitative theories for all the other 
quantities near H C 2- The exact limiting behaviors would 
be useful not only for their own sake, but also for get- 
ting an insight into the behaviors over < B < H C 2 . In 
addition, they would serve as a guide for any detailed 
numerical studies for < B < H C 2 . 

With these observations, I here perform revised cal- 
culations for the Maki parameters and the pair poten- 
tial near H C 2- Besides correcting Eilenberger's results 
on K2 for the spherical Fermi surface, I also perform two- 
dimensional calculations of the quantities for an isotropic 
(i.e. cylindrical or circle) Fermi surface. Thereby clarified 
will be a rather large dependence of K\ (T) and K2 (T) on 
detailed Fermi-surface structures. Indeed, even the em- 



pirical inequality K2 > «i > kgl will be shown violated 
in some cases in two dimensions, even without the spin 
paramagnetism— 

The starting point adopted for these purposes is the 
quasiclassical Eilenberger equations— As emphasized by 
Eilenberger— and also by Serene and Rainerji 6 . the qua- 
siclassical equations have an advantage over Gor'kov 
equations^ 7 - that they are easier to solve due to the ab- 
sence of an irrelevant energy variable. They have a 
rigorous microscopic foundation and hence form a firm 
basis for any quantitative description of superconduct- 
ing/superfluid Fermi liquids. 16 Thus, it seems somewhat 
surprising that few calculations on «2 have been per- 
formed based on the Eilenberger equations— 

This paper is organized as follows. Section ITT1 provides 
the formulation for the s-wave pairing with an isotropic 
Fermi surface and s-wave impurity scattering, deferring 
p-wave impurity scattering to Appendix^] The main re- 
sults are given in Sees. Ill El and ITTFl and the differences 
from Eilenberger's calculation are explained in Sec. Ill HI 
Section IiTlI presents numerical results. Section ]IV\ sum- 
marizes the paper, with possible extensions to include 
realistic Fermi surfaces from first-principles calculations 
and/or anisotropic pairings. Appendix lAl derives an an- 
alytic expression for the magnetization. 



II. FORMULATION 

A. Eilenberger equations 

I consider the s-wave pairing with an isotropic Fermi 
surface and s-wave impurity scattering in an external 
magnetic field H || z. The vector potential in the bulk 
can be written aa 19 : 20 : 21 : 22 : 23 : 24 



A(r) = Bxy+ A(r) 



(1) 



where B is the average flux density produced jointly by 
the external current outside the sample and the supercur- 
rent inside it, and A expresses the spatially varying part 



2 



of the magnetic field satisfying J Vx Adr = 0. I adopt 
the units where the energy, the length, and the magnetic 
field are measured by the zero-temperature energy gap 
A(0) at H = 0, the coherence length £o = hyp / A(0) with 
vf the Fermi velocity, and Bq = 4>q/2tt^q with (f>o = hc/2e 
the flux quantum, respectively. 1 also put h — = 1 
and use the gauge where V ■ A = 0. The Eilenberger 
equations^ now read 



A(r) In 



M 

2r 



■■(V-iA) 



A +^r).'/- (2a.) 



-2vrT 



n=0 



(/(£„, k F ,r)) 



A(r) 



V 2 A(r)= -4,2ttT V(v ff ( £ „,k F ,r)> 



(2b) 



(2c) 



Here £ n = (2n+ l)7rT is the Matsubara frequency, r is the 
relaxation time in the second-Born approximation, (• • • ) 
denotes the Fermi-surface average satisfying (1) = 1, A(r) 
is the pair potential, and the unit vector v = k specifies a 
point on the isotropic Fermi surface. The quasiclassical 
Green's functions / and g are connected by g=(l—ff if ) 1 ^ 2 
with f^(e n , kp, r) = /*(e n , — kp, r), and the dimensionless 
parameter kq is defined by 



(3) 



where H c (0) = ^4ttN(0)A(0) is the thermodynamic crit- 
ical field at T = with N(0) the density of states per spin 
and per unit volume. Equations l|2ajl - l(2"c)l are to be solved 
sclf-consistently for a fixed B. Finally, the missing con- 
nection between H and B is obtained by applying the 
Doria-Gubernatis-Rainer scaling 25 to Eilenberger's free- 
energy functional^ The details are given in Appendix 
A. The final result is given by 



*T ^ f / /ty.(V-,A)/-/v-(V- 



iA)/t 



where V is the volume of the system. 



(4) 



B. Expansion near H C 2 

Near H C 2, the coupled equations J5J and J1J are ex- 
panded in terms of A(r) as follows. First, let us rewrite 2 ^ 



(V - iA) 



Bsin9 r 
2^2 



ltp (a+A)-e iip (^ + A*) 



(5) 



where (6, <p) are the polar angles of v, and the quantities 
a, a' , and A are defined by 



f 



d .d 

1— + Bx 



/2B \dx dy 
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— — — + 1— + Bx 



2B V ® x ®y 



7 . A x + iA y 

A = —i , 

V2B 



(6a) 



(6b) 



with [a, = 1. The operators (a, a^) are the same as 
(a_,a+) introduced by Helfand and Werthamerf 5 - and 
(F-,F+) by Eilenberger >i£ 
I then expand /, g, and A up to the third order in A(r) 

as 



/ = / (D +/ (3) 

ff = 1 _l / ( l )t / (D . 

A=A<$ 



(7) 



Substituting Eqs. J5J and Q into Eq. lj2*a|) and collecting 
terms of the same orders, we obtain 

[e n + /3 (e-^a - e^o*) ] = A + , (8a) 



[e n + (3 (e-^a - e*<f) ] / (3) = 

/ (1)t / (1) /a , </^>\ , (1) 
2 ^ 2r J 4t 1 

-/3 (e-^i^ - e ^l^*)/W , 



with 



Also, Eq. I|2bjl is transformed into 



£ sin 9 
2V2 



A(r)ln^ = -27rT^ 



T 



n=0 



(/ (D) + (/ 0)) _ 



(8b) 



(9) 



(10) 



The Maxwell equation l|2"cj) is given in the leading order 
by 



^(/Wt/We^sin^), (11) 



whereas Eq. ||2J becomes 
ttT 



H - B 



2V2B 



-e itp a^)f (1) sme). (12) 
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C. Two-dimensional calculations 

To investigate the dependence of the Maki parameters 
on the Fermi-surface structure, calculations will also be 
performed for a two-dimensional system with an isotropic 
Fermi surface placed in the xy plane perpendicular to B. 
The analytic expressions for this case can be obtained 
from those of the three dimensions by simply putting 
sin0— > 1 and omitting the integrations over 9. 



Equation l|15|l has a wide range of applicability over 
B > 0.1H c2 both near T c and in the dirty limit. How- 
ever, the region in the clean limit shrinks as T — > to 
disappear eventually. It should also be noted that higher 
Landau levels oiN = even become relevant for anisotropic 
pairings and/or anisotropic Fermi surfaces at low temper- 
atures. 

Substituting Eqs. (|13bf) . ljT3cJl . and JTHJ) into Eq. © 
and using the orthogonality of e lmv and tpNq, we realize 
that f^ N can be written as 



D. Transformation into algebraic equations 

Equations JSJl and (|TUj> - ((T2)l are solved with the 
Landau-level expansion (LLX) method 24 by expanding 
A, (y = 1,3), and in terms of periodic basis 
functions of the flux-line lattice as 



A(r) = W J2 A n ^jv q (r) . 



(13a) 



N=0 



JmN — °rnN^ J N 



(16) 



Equations (JSJ and (|1(J[1 - H12I) are thereby transformed into 
algebraic equations for f( 3 \ Aq, Aj^ , and H—B as 



(i) 

N' 



?N0 



N 1 



(ft 

2r 



(17a) 



/(")(£„, k F) r) = W 



oo oo 



m=—oc N—0 



<••• V'Arq(r), 
(13b) 

(13c) 



Here N denotes the Landau level, q is an arbitrary chosen 
magnetic Bloch vector characterizing the broken transla- 
tional symmetry of the flux-line lattice and specifying 
the core locations, and K is a reciprocal lattice vector 
of the magnetic Brillouin zone. See Ref. 24 for the ex- 
plicit expressions of the basis functions ^>jvq(r); they are 
essentially equivalent to Eilenberger's %pN(r\ro)H and re- 
duces for iV = to Abrikosov's solution for the Ginzburg- 
Landau equations near H C 2*& It only suffices to know the 
properties: 



(lpNq\lp 



N'q, 



= 5 



NN' 



(14a) 



(14b) 



lq 



On the other hand, the expansion in K in Eq. I|13cfl was 
introduced by Brandt*^ for solving the Ginzburg-Landau 
equations over < B < H C 2 . This expansion enables us 
to integrate the Maxwell equation appropriately so that 
J V x A dr = is satisfied automatically. 

A simplification results in Eq. 1)1 3JI for the s-wave pair- 
ing with an isotropic Fermi surface near H C 2- Indeed, 
A(r) can be described excellently with only the lowest 
Landau level asS4 



(15) 



A(r) = %/T/A Voq(r). 



N' 



(ft 



(3)\ 



2r 



r(3) 
'N 



'N ) 



n=Q 



A 



(2) _ 27tTAq ^ 



K 



EE 

n=0 N 



J%\e n )I N+1N {K) 







(17b) 



In | = -27rrf (</«) + (f^)Al 1) , (17c) 



, (17d) 



(17e) 



n=0 AT 

Here the matrix M. is defined by 

Mnn> = Snn>£u + 6 n ,n'-i0VN+1 - Sn.n'+iP^ , 

(18) 

and Jn's are given by 

jff = t^y/W+l { fW f^ ) , (19a) 



J 



(3) 
N 



2 V 



B 



(-1)" J 



N' r (4) 7(1) 
NN'N+N'O J N' JN+N 




^ T Z_^\ L ) 1 NN'N'N \JN' J N' IJN ' \ lmJ ) 



'N 



N' 



JL 

An 



E [^N+IN (^)^k In+1 



K^O 



-Inn -i (K) A^ * fffi_ J 



(19c) 
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with 2 ! 



l( NlN2N 3 m = V I V^qV^qVWaqV^qdr , ( 2 0a) 



I NlN . 2 (K) = y^ iq ^ 2q e- lK - r dr. (20b) 



Equation (|17a() tells us that fy' is real; this fact has 
been used in wr iting down Eqs. llTEI) - ((T7el) . As for /^ 3) , 
numerical calculations show that I^'s appearing in Eq. 
(1 19 bp are all real for the relevant hexagonal lattice, with 



/?a ^41 -1.16. 



(21) 



Also, Jjv+iat(K) can be transformed with partial integra- 
tions as 

I N+1N (K) = 7 =L=p / ( a t^ q )> Wq e- jK - r rfr 
V Jv + 1 J 

JjVjV_l(K) + 7== = /jyjv(K) 



Ky-iK x 



^2B{N+l) 



A' 



V2B(iV+l)^ 



(22) 



Since -Zjvat(K) is real, J'' 4 ) is also real from Eq. (|17d(l . 
and so is fffi. 

It is desirable for a later purpose to express J^ A > in 
terms of 1^ rather than 7jv + iat(K). This can be per- 
formed by first substituting Eq. I|17d|) into Eq. (|19cl) , and 
then using Eq. (f22J an d the identity J2k^o e * K '■ r ~ r ' ) = 
VS(r—r') — l. The result is given by 



(A) _ irT(3 



N 



]T J N> [VN+lfSilPKN' ~ 1) 



u n'=0 AT' 



AT/^Xjv-ijv'-I)], (23) 

where /^4i = /jvilfcn)' J$ = J N>( £ n>) is given by Eq. 
(|19a|l . and Iaw is an average of l' 4 ^ defined by 



Inn' 



1 



TV JV' 



{N + \){N' + l) 



E E 

(24) 



JVi=0 JV,'=0 



E. Solutions 

We are now ready to solve Eqs. (|17afl and l|17b|l . To 
this end, let us define 

\N+N' f>N 



(25) 



where the second equality originates from Mnn 1 — 
{-l) N+N ' M N 'N- Then Eq. ((T7a|) is transformed into 



Solving Eq. (|26|l self-consistently for (/q 1 ^) and substi- 
tuting the result into Eq. H26|) . we obtain 



JV 



l-«/2r 



(27) 



The denominator in Eq. I|27|l corresponds to the so-called 
"vertex correction." Equation (|17b[) for fffi may be han- 
dled similarly. Using the symmetry Kq = (— 1) N K^, we 
thereby arrive at the expression for the relevant quantity 
(fk 3) ) in Eq. (JlTcJ) as 



(28) 



JV 



The quantity Ejvv -1 )^//^ 4^) in E Q- © ma y be 
transformed further by using Eq. H23[) and (|19af) as 

^iH/I^oj^M) 

V 

3 £ J^M 2ttT£ E ^(^)(W-1) • (29) 



JV 



JV 



=0 JV' 



We next consider the self-consistency equation (|17cf) 
for the pair potential. Here, fjp (B) is expanded in terms 
of the distance H c2 — B from H c2 as 

fV(B) = f ( N ] (H c2 ) - f^'(H c2 )(H c2 -H+H-B) , (30) 

whereas the higher-order term is evaluated at H c2 . 

To find an explicit expression for fff' in Eq. H30(l . let us 
differentiate Eq. (|17af> with respect to B: 



EM f - A ^° 



(1)'\ 



JV' 



2r 



(i) 



(31) 



This equation can be solved in the same way as Eq. I|17a|) 
to yield 



</n— e^ 



(2) 
JV ' 



(32) 



N 



where jffi is defined by Eq. (|19a|l . Let us substitute 
Eqs. and (E3 into Eq. JTZ3), replace H-B by 

the right-hand side of Eq. (|17ef> . and regard H c2 — H as 
second order. Collecting first-order terms, we obtain the 
equation to fix the upper critical field H c2 as 



In 



To 
T 



(33) 



The third-order terms determine the pair potential Ao 
and the magnetization H—B as a function of H c2 —H as 



Jn — a jv 



2r 



(26) 



A 2 - 



H c2 — H 



4Sa/S% - Sa/SI s 2 ' 



(34a) 



5 



H -B 



H c2 ~H 



H c2 — H 



where £2, £4, and Sa are defined by 



71=0 N=0 



(34b) 



(35a) 



X>iV = det 



They satisfy 



in VTP 



-VN^IP in 



V N -t = i n V N + N[i 2 V N+1 , 



(37b) 
(38a) 



S 4 



(3), 



n=0 N=0 



00 00 



00 



00 



s A = (2ttt) 2 EEEE jgHen^fayrNN' , 



E 

n=0 n'=0 iV=0 iV'=0 



(35c) 

with /j^ , Jjv, and 1nn> given by Eqs. (f2*7j) . and 
()24p. respectively. All the quantities in Eq. I|35|l are to 
be evaluated at H c2 . The latter equality in Eq. I|34b[) 



defines the Maki parameter k 2 with 0& = L 



(4) 
0000 



= 1.16. 



Equation H34|) forms the main result of the paper, 
which is not only exact but also convenient for numeri- 
cal calculations. An extension to include p-wave impurity 
scattering is carried out in Appendix B, where it is shown 

that Eq. I|34|) is still valid with the replacements of 



and jffl by Eqs. JES) and ljB9|) . respectively. 

Sometimes it is physically more meaningful to express 
Ao as a function of B instead of H , because B is the real 
average field inside the bulk directly relevant to the spa- 
tial profile of the pair potential. It is obtained, without 
the replacement of H— B mentioned above Eq. (J33J), as 



A 2 
^0 



Hrl — B 



4S4/SZ - Sa/SI 



1 S 2 



(36) 



f>N' 



F. Calculation of K 

The key quantity in Eq. is K§' defined by Eq. 
(J2SJ), as may be seen from Eqs. l|2"7)l. ifT^l) . and An 

efficient algorithm to calculate them is obtained as fol- 
lows. 

Let us define T>m (Pn) for N = 0, 1, 2, • • • as the de- 
terminant of the submatrix obtained by removing (re- 
taining) the first ./V rows and columns of the tridiagonal 
matrix Ai of Eq. I|18fl , namely, 



T>n = det 







i n VN+l/3 

■y/W+\P in y/N + 2(3 







-2/3 



(37a) 



V N+1 = i n V N + N(3 2 V N _x , 



(38b) 



as shown by expanding Eqs. Ij37a(l and l|37bjl with respect 
to the first and last row, respectively^ Then for 
N' <N is obtained by also using standard techniques to 
solve linear equations21 as 



K 



N' 
N ■ 



N—N' 




N\ V N+1 V 



N' 



(39) 



This algorithm can put into a more convenient form in 
terms of 



They satisfy 



Ti-N = inT^N+i/I^N > 



TZn — i n T>N^i/T>N 



TZn-i 



l + Nb 2 Tl N ' 



f 



1 + Nb 2 K 



N 



with IZi = 1 and 



b = PI in ■ 



Then K% and K%' for N'<N are obtained by 



K° = K /in 



K% = {K N /1l N )K^Zi , 



N-l 



Kg = VNbKNKfiU 



(40a) 



(40b) 



(41a) 



(41b) 



(42) 



(43a) 



(43b) 



(43c) 



Numerical calculations of IZn may be carried out by 
starting from 72-jv cllt = 1 for an appropriately chosen large 
2V cu t and using Eq. (|41a|) to decrease N. One can check 
the convergence by increasing N cut . It turns out that 



6 



N cu t — 1 is sufficient both near T c and in the dirty limit, 
thereby reproducing the analytic results by Gor'ko\ii£ 
and Caroli, Cyrot, and de Gennes^ respectively. In 
contrast, -/V cut > 1000 is required in the clean limit at 
low temperatures. 

Noting that Eq. with Eq. J2ZJ) should be equiv- 
alent to the equation for H c2 obtained by Helfand and 
Werthamer^ we get an alternative expression for as 



e' x2/2 dx . (44) 



IT J Q £? + X 2 f3 2 



The equivalence between Eqs. (|43a|l and (|44[> can also be 
checked numerically. 



G. Expression of kgl 

Near T c where (3 <C e n holds, we may choose iV cu t = I 
in Eqs. I|41|) - (|43l) and expand the resulting expressions 
with respect to f3/e n - This yields Kq w f/e„ — /3 2 /e^ and 
Ki ss /3/e 2 , so that Eq. l(77|) can be approximated by 



/o (1) « l/e„-(/3 2 e, i + (/3 2 )/2r)/ e 2 £ l and ft 1 ' « /3/ E „e„. 
Using these results in Eq. (|35[1 and retaining only terms 
of the leading order in /?, we obtain 



(i) 



2d(nT c ) 2 ^ (2n + l) 2 (2n+l + l/27rrT c ) 



5 4 



7C(3) 
8(^T C )< 



-5.4 = S 2 0A 



(45a) 



(45b) 



(45c) 



where d = 2, 3 is the dimension of the system. Substi- 
tuting Eq. 1451) into Eq. (|34b|) . we find the expression of 
kgl = k 2 (T c ) as 



d7rT cV /7C(3)/2 



^[(2n + l) 2 (2n+l + l/2 7 rrT c )]- 1 



k . 



(46) 



This expression enables us to eliminate kq in favor of 
kgl- 

The case with p-wave impurity scattering may be 
treated similarly by using Eqs. (SU and ||E9j| for ffl 

(3) 

and Jjy , respectively. The resulting kql is given by Eq. 
(|46[1 with a replacement of r by the transport lifetime re- 
defined through 



f 

Tu- 



rn agreement with Eilenberger 



I 

n 

12 



(47) 



H. Eilenberger's results 

I now clarify the connection with Eilenberger's well- 
known results— They are obtained by extracting from 
lff N N N of Eq. l|2Ua|) a part which may be expressed 
in terms of /3 A = 1.16 of Eq. JSTJ. 

To see this, let us start from an alternative expression 



of/ 



(4) 

N!N 2 N3N4 



for Ni+N 2 =N 3 +Ni: 



J 



(4) 

N 1 N 2 N 3 N 4 



2 



\lpN a O a (0) 



N, 



x(NiN 2 \Ni 



-N 2 -N a N a )(N 3 N 4 \N 1 +N 2 -N a N a ) . 

(48) 

Here ipN a Oa(0) and (NiN^N^N^ are the quasiparticle 
wave function and the overlap integral defined by Eqs. 
(3.12) and (3.23) of Ref. El respectively. This identity 
can be proved by using Eqs. (3.22) and (3.24) of Ref. 

IH and noting ^^(r) = ip < Nq{'^ r ) which both denote the 
present ?/>jv q (r). If we retain only terms of N a — in 
Eq. and use ^ ELi IV'ooa(0)| 2 =j3x, we obtain an 
approximate expression for Eq. (|48|l as 



r (4) 

L N 1 N 2 N 3 N 4 



/9A<JVlJV2|JVl+JV20) (N 3 N 4 \N r 

(Ni+N 2 )\ 



-N 2 0) 



2^+^ y/N^.N^.NslNi 



I 



(4E) 



l N 1 N 2 N 3 N i ■ 

(49) 

Now, Eilenberger's approximation for k 2 is given by Eqs. 
(2.7), (2.8), and (6.5) of his paper— and corresponds to 



k 2 « (K 2 Sr } /2S 2 l3 K ) 

j( E ) 



(50) 



1/2 = (E) 



in Eq. (|34b|) . where S4 is obtained from Eq. H35b|) by 

replacing J*- 4 ) by /( 4E ) . Indeed, this procedure yields nu- 
merical agreements with his results. As seen below in Sec. 
however, this approximation is not correct quanti- 

■1%. 



imBl 

tatively far beyond his estimation 

It should also be noted that Eilenberger's definition of 
k 2 by Eq. J3U|) is different from Maki's through Eq. (|34bl) 
with respect to 



77 = S A /S 2 (3 A 



(51) 



I resume Maki's definition through Eq. (|34b|) where k 2 
has a one-to-one correspondence with the initial slope of 
the magnetization. This is certainly more preferable than 
expressing the slope with two parameters, k 2 and 77. 

I finally comment on Eilenberger's analytic expression 
for 77. In addition to the approximation 149|) . it was ob- 
tained by integrating the Maxwell equation with a re- 
moval of a common operator; see the argument above 
Eq. (6.4). 12 However, this procedure may bring an erro- 
neous constant. Indeed, his B (r) of Eq. (6.4) does not 
satisfy the required condition J B (r) dr = 0. Thus, his 
expression for 77 is incorrect in the two respects and can- 
not be obtained by adopting the approximation 149() in 
Eq. JSTJ. 
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FIG. 1: Temperature dependence of k^/kgl for several values 
of £&/kr in the extreme type-II case acgl = 50. (a) l tT /l = 1.0; 

(b) kr/l = 2.0. 



III. NUMERICAL RESULTS 



A. Numerical procedures 



I have adopted the same parameters as those of 
EilcnbergeriiS 



&/ltT = l/2irT cTtI , kr/l = rtr/r , (52) 



FIG. 2: Temperature dependence of r\ defined by Eq. I5H for 
several values of £&/ki with k T /l = 1.0. 



B. Results for K2 

Figure 1 shows K2/kql as a function of T/T c for dif- 
ferent impurity concentrations. The upper one is for 
ltr/1 = 1-0, i.e., the case without p-wave impurity scat- 
tering, whereas the lower one is for l tI /l — 2.0. They 
are calculated in an extreme type-II case of kql = 50, 
so that they directly correspond to Eilenberger's re- 
sults for l tr /l = 1 and 2, respectively^ These curves 
show qualitatively the same behaviors as those of Eilen- 
berger's, including the divergence in the clean limit for 
T 0, as predicted by Maki and Tsuzuki£ Except the 
curves in the dirty limit, however, marked quantitative 
differences are seen. For example, K2(T = 0)/kgl for 
(^E/itr, hr/l) = (1-0, 1.0) is 1.40 from the present calcula- 
tion, whereas it is 1.50 from Eilenberger's. Thus, we real- 
ize that Eilenberger's approximation (|50|l yields quanti- 
tative errors of < 20% for the deviation K2/KGL— 1- Com- 
paring the two figures, we observe the followings: (i) The 
results in the dirty limit are the same between lt T /l = 1 
and 2; (ii) p-wave scattering has a general tendency to 
lower the values of K2, and also produces a non-monotonic 
behavior in K2/kql as a function of ^/hr- 

Figure 2 displays rj defined by Eq. (|51|l as a function of 
T/T c for £ E /7tr = 0.05-50.0 and l tT /l = 1.0. This quantity 



to express different impurity concentrations. Numeri- 
cal calculations of Eqs. I|33(l - (|35[1 have been performed 
for each set of parameters by restricting every summa- 
tion over the Matsubara frequencies for those satisfying 
£ n < s c . Choosing e c — 50^ 100 has been sufficient to ob- 
tain an accuracy of ~0.01% for K2- On the other hand, 
summations over the Landau levels have been truncated 
at N = N cut where I put TZjy cu% = 1 in the calculation of 
K^j ; see Sec. Ill Fl for the details. Enough convergence 
has been obtained by choosing iV cut = 4, 40, 100, 200, 500, 
and 2000 for ^ E /Z tr = 50, 1.0, 0.5, 0.1, and 0.05, respec- 
tively. Finally, integrations over 9 have been performed 
by Simpson's formula with N cnt +1 integration points for 
O<0<tt/2. 
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FIG. 3: Temperature dependence of K2/kgl for kgl = 1 with 
kr/l = 1.0. 
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FIG. 4: Temperature dependence of k^/kgl for an isotropic 
two-dimensional system in the extreme type-II case kgl = 50. 
(a) l u /l = 1.0; (b) kr/l = 2.0. 



becomes relevant for small values of kgl at low temper- 
atures, as may be realized by Eq. (|34b|) . The curves also 
deviate substantially from Eilenberger's results. For ex- 
ample, r) for ^E/^tr = 0.25 at T — is 1.34 from the present 
calculation, whereas it is ~ 1.11 from Eilenberger's. Gen- 
erally, the values are larger than those of Eilenberger's. 
This fact implies that Ki (T) for kgl ~ 1 becomes smaller 
than the evaluation of Eilenberger. 

To see the dependence of K2/KGL on kql explicitly, I 
have performed a calculation of near the type-I-type- 
II boundary of kql = 1.0. Figure 3 plots the results for 
£ E /Z tr = 0.0-50.0 and l tT /l = l. Compared with Fig. 1(a), 
we observe that each curve is slightly shifted downward. 
However, the changes are surprisingly small, considering 
the closeness to the type-I-type-II boundary. We thus 
realize that the factor Sa/ S% = ??/3a in Eq. (|34b|) can be 
neglected practically for kgl>5, as already observed by 
Eilenberger iiS 

The above calculations are performed for an idealized 
spherical Fermi surface. However, real superconductors 
are often characterized by complicated Fermi surfaces. 
To see the dependence of K2/KGL on Fermi-surface struc- 
tures, I have performed an isotropic two-dimensional cal- 
culation described in Sec. Ill CI Figure 4 shows the re- 
sults, where the parameters are the same as those in Fig. 
1. The curves for ^E/^tr = 50 are almost the same as those 
in Fig. 1. Thus, in the dirty limit, we have a universal 
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FIG. 5: Temperature dependence of ki/kgl for several values 
of £ E /Ztr with l tr /l = 1.0. (a) d = 3 ; (b) d=2. 



curve which depends neither on detailed Fermi-surface 
structures nor fine features of the impurity scattering. 
As the system becomes cleaner, however, differences due 
to the two factors emerge eventually. In fact, we observe 
that each curve for £e/^i-<1-0 in Fig. 4 deviates far less 
from 1 than the corresponding one in Fig. 1, and the 
temperature dependence is also weaker. Another point 
to be mentioned is that, even for ^/hv = 0.05, we see 
no trace of divergence as T — > 0. Indeed, a closer ex- 
amination of the analytic results by Maki and Tsuzuk: 8 
and Eilenberger— enables us to realize that it is the re- 
gion 9 ~ in three dimensions which is responsible for 
the divergence of k>2- Thus, we may conclude that K2 in 
two dimensions remains finite even in the clean limit as 
T — ► 0. In general, Ki will remain finite if the relevant 
Fermi surface does not close along the direction of the 
magnetic field. 



C. Results for k\ 

The Maki parameter k\ is defined by 2 



K\ = H, 



c2 



2H r 



(53) 



where H c = H C {T) is the thermodynamic critical field. 
The preceding results for suggest that ki(T)/«gl m ay 
also exhibit considerable dependence on detailed Fermi- 
surface structures. 
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Figure 5 compares Ki(T)/kgl between two and three 
dimensions for ? tr /Z = 1.0. The curves for fa/hi = 50 show 
almost the same behavior. As £e/^i- becomes smaller, 
however, the two cases display a marked difference. In- 
deed, ki(T)/kgl is seen to increase (decrease) in three 
(two) dimensions as £&/lt T — *0. 

Figure 6 shows curves of k± (T) / kql in two and three 
dimensions for lt T /l — 2.0. Again the p-wave impurity 
scattering is seen to lower the value of ki/kgl, and also 
introduces non-monotonicity in ki/kql as a function of 
£,E/kr- Especially in two dimensions for ^/hr = 0.1-1.0, 
ki/kgl becomes smaller than 1 over finite temperature 
ranges, i.e., the emperical inequality k,i > K\ > kgl is not 
satisfied here, even without spin paramagnetism^ 

A substantial dependence of H C 2 on Fermi-surface 
structures may be realized more clearly by looking at 
the temperature dependence of the reduced critical field 
introduced by Helfand and Werthamer^ 



h*(t) 



H c2 (t) 



-dH c2 (t)/dt\ t=1 



(54) 



where t= T/T c . Figure 7 compares h*(t) between two 
and three dimensions for both the clean and dirty lim- 
its. The curves coincide in the dirty limit, whereas those 
in the clean limit show a marked quantitative difference. 
We also observe that h* (t) in two dimensions is a rather 
sensitive function of purity. A considerable reduction of 
h* d=2 (t) in the pure limit from h* d= ^(t) may be attributed 



FIG. 7: Temperature dependence of the reduced critical field 
h*{t) for the dirty limit of d = 2, 3 (mid curve) and the clean 
limit of d = 2 (lower curve) and d~3 (upper curve). 



to the pair breaking by supercurrent. This effect is more 
effective in two dimensions. Indeed, a point on the cylin- 
drical Fermi surface is equivalent to a point on the equa- 
tor of the spherical Fermi surface perpendicular to H 
where the pair breaking is most effective. This fact can be 
seen clearly in the polar-angle dependence of the density 
of states calculated by Brandt, Pesch, and Tewordt. 29 
Put it another way, if the relevant Fermi surface does 
not have a closed orbit perpendicular to H, the corre- 
sponding h*(t) in the clean limit will be enhanced over 
the prediction for the spherical Fermi surface. 

A considerable reduction of h*(t) or K±(t) in the 
presence of spin paramagnetism was established by 
Werthamer, Helfand, and Hohenberg, 6 and also by 
Maki. 11 The present results indicate unambiguously that 
the Fermi-surface structure is also an important fac- 
tor for h*(t) in clean systems, as already noticed by 
Helfand and Werthamer, 5 Hohenberg and Werthamer, 30 
and Werthamer and McMillan^ 



D. Results for the pair potential 

A quantity of fundamental importance is the coeffi- 
cient A , which is equal to the spatial average -J (| A(r)| 2 ) 
of the pair potential and relevant to all thermodynamic 
and transport properties near H C 2- It is physically more 
meaningful to express it as a function of the real average 
field B in the bulk instead of H . Equation shows 
that Aq(B) is proportional to (H C 2 — B) 1 / 2 near H C 2- I 
here express this A by using the energy gap A(T) at 
5 = as 



A (B,T) = c (T)(l- J B/i7 c2 ) 1 / 2 A(T) 



(55) 



Then the coefficient c(T) should be of the order of 1. 

Figure 8, calculated for the spherical Fermi surface, 
displays temperature dependence of c(T) in an extreme 
type-II case of kgl = 50 for (a) l tT /l = 1.0 and (b) 
1^/1 = 2.0. Thus c(T)~l, as expected, having the same 
value 0.929 at T c . Differences among different ^/hi grow 
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in the extreme type-II case kgl = 50 as a function of T/T c for 
several values of (,E/hr- (a) hi/l = 1.0 ; (b) = 2.0. 



at lower temperatures, and c(T) for £E/7tr < 0.1 drops 
rapidly near T = 0. Indeed, c(T) in the clean limit for 
three dimensions is expected to reach as T — ► 0, cor- 
responding to the divergence of ^2. This also implies 
that the expansion in A(r) near H C 2 is no longer valid in 
this limits* The curves in the dirty limit are the same 
between l tI /l = 1.0 and Z tl /Z = 2.0. For Wtr<l-0, how- 
ever, each curve for lt T /l = 2.0 at low temperatures has 
larger values than the corresponding one for lt T /l = 1.0. 
Thus, finite p-wave scattering in clean systems tends to 
increase c(T). 

The coefficient c(T) also increases mildly as kql be- 
comes smaller, as realized by comparing Fig. 9 for kql = 1 
with Fig. 8(a) for k G l = 50. 

Figure 10 plots results of the two-dimensional calcu- 
lations performed with the same parameters as those 
in Fig. 8. The curves for the dirty limit are the same 
between two and three dimensions. As the system be- 
comes cleaner, however, the coefficient c(T) for two di- 
mensions becomes larger than the corresponding one for 
three dimensions. Thus, for clean systems, we observe 
once again a considerable dependence of the coefficient 
c(T) on Fermi-surface structures. 



IV. SUMMARY 

This paper has presented revised calculations of the 
Maki parameters «i and Ki as well as the spatial average 
(|A(r)| 2 ) near H C 2- Eilenberger's results for K2 have been 
corrected appropriately, as described in Sec. Ill HI The 
analytic expressions derived in Sees. Ill El and III Fl have 
been useful to carry out efficient calculations for both two 
and three dimensions with isotropic Fermi surfaces and 
arbitrary impurity concentrations. Thereby found are 
large quantitative differences of the parameters between 
two and three dimensions (except in the dirty limit where 
there are no differences between the two cases). For ex- 
ample, no trace of divergence in K2(T — > 0) is found for 
the clean limit in two dimensions. 

The present results clearly indicates the necessity of 
considering detailed Fermi-surface structures from first- 
principle calculations for an quantitative understanding 
of the Maki parameters in clean superconductors. This 
was already recognized by Helfand and Wcrthamer^ by 
Hohcnbcrg and Werthamer^ and also by Werthamcr 
and McMillan 31 when their strong-coupling calculation 
could not explain a large deviation of k,\ / kgl observed 
in pure niobium22i2£ and vanadiumM from the theoreti- 
cal prediction of Helfand and Werthamer. 5 Efforts have 
been made along this line to establish a realistic calcula- 
tion of K\, or equivalents, ff c2 , 30 i 35 i 36 i 37 i 38 i 39 i 40 i 41 i 42 i 43 i 44 
However, little progress seems to have been achieved with 
respect to K2- 

The method developed here for k\ and K2 may be 
extended easily to include Fermi-surface structures and 
anisotropic pairings. Some of the necessary modifica- 
tions are: (i) to use the general expansion (|13af> with 
N = even for the pair potential, rather than Eq. I|15D ; (ii) 
to use more convenient basisfunctions than e mlip in Eq. 
(|13bll for describing the kp dependence of /(e n ,kp,r), 
such as the Fermi-surface harmonics of Allen4^S4Li£ 
The corresponding matrix M. in Eqs. I|17a|) and (|17b|) is 
no longer tridiagonal, but may be inverted rather easily 
with present high-speed computers. 
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APPENDIX A: DERIVATION OF EQ. |gj> 

To obtain Eq. QJ, let us start from Eilenberger's free- 
energy functional^ per unit volume with B chosen as an 
independent variable instead of H. It is given in units of 
iV(0)A(0) 2 as 



V 



|A(r)| 2 ln 



|A(r) 



where / is defined by 
1 = A*/ + A/ t + 2f„(.g-l) 



- (/(£„, k F ,r)) 



T 
(Al) 



+(<?-l)v 



/t(V-iA)/-/(V+iA)/t 
2//t 



4r 2r 

(A2) 



The functional derivatives of Eq. IjAlfl with respect to 
/t, A, and A lead to Eqs. |(2a||-|(ScJ, respectively. The 
last term in Eq. (jA2|l is slightly different from the origi- 
nal functional of Eilenberger where gv appears in place 
of (g— l)v. 15 Although it docs not change Eq. (0) at 
all, it is found numerically that the modification is nec- 
essary for F to have its absolute minimum with re- 
spect to A, A, and / = /(A, A) satisfying Eq. l(2"a|) . 
as anticipated by Eilenberger— It should be noted that 
Pesch and Kramer— also adopted Eq. HAlfl as a basis 
for their numerical calculations. More recently, Endres 
and Rainer— performed a numerical calculation of the 
free energy for both an SN contact and a single vor- 
tex based on Eq. IjAlf) . and compared the results with 
those from three free-energy functionals obtained from 
the Luttinger-Ward functional. They found numerical 
agreements among the values from the four different ex- 
pressions. 

Following Doria, Gubernatis, and Rainer^ I now 
rewrite the right-hand side of Eq. (|A1|) in terms of 



r'=r/A, B X = \ 2 B, A A (r') = AA(Ar'), 
A A (r') = A(Ar'), / A ( e „,k F ,r') = /(e n ,k F ,Ar') 



(A3) 



I then differentiate the resulting expression with respect 
to A and put A = 1 . Since the procedure (|A3|) does not 
change the value of F/V, we have J^(F/V)\\ = i =0 from 
the left-hand side. As for the right-hand side, the only 
implicit dependence to be considered is the one from By, 
those from f\, A A , and A\ can be neglected due to the 
stationarity of Eq. (J2J. We thereby obtain 



= 



d{F/V) dB x 



dB x d\ 



^2k 2 B 2 -^ frfr(VxA) 
=i ' J 

nT™ f / /ty.(V-»A)/-/v.(V + zA)/t 



(A4) 

Using the thermodynamic relation d ^ V ^ — KqH in 
the present units, we arrive at Eq. 



APPENDIX B: EXTENSION TO THE CASE 
WITH p-WAVE IMPURITY SCATTERING 

In the presence of p- wave impurity scattering, Eq. Ij2a(l 
is replaced by 



(g) dk-(k'g) v 
2r 2ri 2 



:V - iA) 



= ( A+ ^ + ^r l ) g - 



(Bl) 



where (k'g) = (k'g(e n , k F , r)), for example, and d — 2, 3 
is the dimension of the system. This brings additional 
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terms on the right-hand side of Eqs. i|17a|) and l|17b(l as 
E^W/$ =S m (l 



Jivi- 



N' 



2t 2n 
dsin6»(/ 1 (1) sin6»') 

( 3 )\ ^^ofl/-f( 3 )^ 



(B2) 



' A" 



(/<5 



dcos0(/ o w cos0') 



2t 

(3)o;, 



2n 



dsmf^/^sinf?' 
4^ 



J 



(3) 
AT 



J 



01) 



(B3) 



where J)^ is given by Eq. (Eg, and J^> is defined 
instead of Eq. (|19b|l by 



t(3) = _ 1 V^/ i\N'A4) 7(1) 
w — nZj' ' NN' N+N'O J N' J N+N' 



A" 



(/ (1) ) , dcos0</^cos0') 



X 1 + 



(1) 



2r 



2ti 



4 EC- 1 )"'* 



N' N 



N' 



rising 



n 



./AT 



Sri 



AT' 



-J 



(4) 

NN'N+N 



'-11 J A" JN+N'-l 



x (/«sin^) +4^'iV'-iiv + i (/^/^-iSin^/Vii 



(i)f(i) 



~ tl NN'N' + lN-l \JN< JA r ' + l t,mt7 IJN-1 



;(i) 



r(l) 



(B4) 



Equation i|B2() can be solved in the same way as Eq. 
lUiJ to yield 



2t 2ri 

&1 dsme(fPsme') 



-K 



N " 



4n 



(B5) 



From Eq. ljB5|l . we obtain self-consistent equations for 
(jfW), and (/^sine) as 







" <^o°> " 


(/^cosfl) 




(Kgcoa6) 






_ (K°sm8) _ 



K 



where the matrix K, is defined by 

{Kg) d{Kgcos9) d{Kf sin (9) 



(B6) 



it 



1-- 



2r 



2-1 



4n 



(itgcosd) d{K$cos 2 9) d{K$sm29) 



2t 



2 n 



8n 



(Jrfsin0) d(^sin26>) s d^sin 2 
2r 4ti 4n 



(B7) 



Noting = K%,(e n ,/3) as seen from Eqs. {HJ-fEJ) 
with /3 defined in Eq. @, we immediately realize that 
(.Kg cos 9) = in Eq. (JEfjJl and £ 2j = £ i2 = for j = 1, 3 
in Eq. (IB7P. Thus, Eq. (|B6|) can be solved easily with 



(/^cos 

/ci 



- 0. Substituting the resulting expressions of 



(/ (1) ) and (/fWfl) into Eq. ifESJl. we obtain 



f(i) 

J AT 



(1 - ^ (^sin 2 gQ)g + ^ (jgsin sin g 



(i-£TO)(i-£<*J«k 2 * , >) 



(B8) 



Equation (|B4|) is also simplified into 



t(3) = lv^/ n AT'r(4) ?(!)?(!) A , (jjL_)A 

A 7 — O"^ NN' N+N'O J N' J N+N' I ' 9.7- / 

A" ^ ' 



1 r( 4 ) /ft 1 ) f(!)\ f(!) 

7~ / A ^ 1 NN'N'N \J N' J N' I J N 



+ - 



4r 

N' 

dsm9 



8ri 



EC- 1 )"' 



A" 



f(4) 
NN'N+N 



(1)^(1) 



'-11 /a" /. 



AT+AT'-l 



x (/TWS') + 41v'jv'-iiv+i {WPnU^I 



^- , A r AT'A" + lA r -l WA" ^A" + l blut7 /JN-1 



(1) 



AT+1 



(B9) 



Equation l|B3|l can be treated similarly to obtain (/q 3 ') 
which appears in Eq. I|17c[) . Then a straightforward cal- 
culation leads to the same expression (|2"51) for (f^) with 
/jy and jj^ replaced by Eqs. I|B8I) and (|B9|) . respec- 
tively. Another relevant quantity in Eq. I|17c[l is (f^ ), 
as seen from Eq. (|3U|) . Differentiating Eq. l|B2fl with 
respect to £> and solving the resulting equation self- 
consistently, we also obtain Eq. (|32|) for (/^ ) with 
in Eq. (|19a|l replaced by Eqs. (jB8j). 



It hence follows that Eqs. I|34|l and 135fl remains the 
same in the presence of the p-wave impurity scattering 
with and J^ 3) replaced by Eqs. and CE3>, re- 



spectively. 
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